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The type lib SN 2008ax: the nature of the progenitor 
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ABSTRACT 

A source coincident with the position of the type lib supernova (SN) 2008ax is identified in 
pre-explosion HST WFPC2 observations in three optical filters. We identify and constrain two 
possible progenitor systems: (1) a single massive star that lost most of its hydrogen envelope 
through radiatively driven mass loss processes, prior to exploding as a helium-rich Wolf-Rayet 
star with a residual hydrogen envelope, and (2) an interacting binary in a low mass cluster 
producing a stripped progenitor. Late time, high resolution observations along with detailed 
modelling of the SN will be required to reveal the true nature of this progenitor star. 

Key words: supernovae: individual: SN 2008ax - galaxies: individual: NGC 4490 



1 INTRODUCTION 



Although the progenitors of several type II supemova e (SNe) have 
been i dentified in pr e-explosion observations (e.g., ISmartt et al.l 
|2004 iLi et al.l |2006|) , the stripped stars thoug ht to explode as 
type Ib/c SNe have so far e l uded discovery ( e.g., Maund & Smarttl 
2004 iMaund et alj 120051 : iGal-Yam et al] 120051 : ICrockett et ail 



20071 l2008h . Thestars exploding as type lib SNe (those events 



which transition from type II to type lb) are believed to be inter- 
mediate of the hydrogen-rich (H-rich) supergiant type II and the 
H-free type lb precursors. The progenitor of the nearby type lib 
SN 1993 J was shown to be such a star; a K0 supergiant stripped 
of most of its hydrogen envelope by a massive binary companion 
dAldering et alj'l994lMaund et alj2004h . 

In this letter we report the identification and characterisa- 
tion of the progenitor of the type lib SN 2008ax, in the galax 



NGC 4490. SN 2008ax was discovered by iMostardi et alj J20oi 
on 2008 Mar 3.45 at a 200 o = 12 h 30 m 40. s 80, <5 2 ooo = 41°38'14."5, 
just 6 hrs after its non-detection in an unfiltered observation made 



by R. Arbour with a limiting magnitude of 18.5 jNakand 120080 . 
SN 2008ax was spect roscopically classified as a type lib SN by 
Chorn ock et al. I J2008h . We have obtained extensive and high qual- 
ity spectro-photometric monitoring of t his SN, the analysis of 
which is presented in a companion paper JPastorello et alj2 008). 



2 OBSERVATIONS AND DATA REDUCTION 

The pre- and post-explosion observations analysed in this study are 
summarised in Table [T] Archival Hubble Space Telescope (HST) 
images of the site of SN 2008ax prior to explosion were recov- 
ered from the STScI archivj] and calibrated via the on-the-fly- 
recalibration (OTFR) pipeline. These observations were taken us- 
ing the Wide Field Planetary Camera 2 (WFPC2) through four op- 
tical filters on three different epochs. The site of SN 2008ax was 
located on the WF3 chip of the F450W and F814W frames, the 
WF2 chip of the F606W image and the PC chip of the F300W im- 
age. The WF chips have pixel scales of 0.1" pix~' while the PC 
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Table 1. Pre-and post-explosion observations of the site of SN 20()8ax. 



Date 


Telescope/Instrument 


Filter 


Exp. Time 








(s) 




Pre-explosion ima; 


;es 




2001 Jul 2 


HST/WFPC2 


F300W 


600 


2001 Nov 13 


HST/WFPC2 


F450W 


460 


1994 Dec 3 


HST/WFPC2 


F606W 


160 


2001 Nov 13 


HST/WPFC2 


F814W 


460 




Post-explosion images 




2008 April 3 


Gemini-N/Altair+NIRI 


K 


1200 



chip pixel scale is 0.046" pix~'. Each of the WFPC2 observations 
were taken as two exposures which were combined in order to re- 
move cosmic rays. Pixel offsets noticed between the F606W ex- 
posures were taken into account during the combination process. 
Finally the combined images where corrected for geometric dis- 
tortion. These calibrated images were used to perform the image 
alignment described in §3. PSF photo metry was carrie d out using 
the hstphot package (version 1.1.7b) (Dolphin 2000a). Option 10 
was chosen, which turned on local sky determination and turned off 
aperture corrections since there were no good aperture stars. (In this 
case hstphot applies default filter-dependent aperture corrections). 
The F606W exposures were input as individual frames along with 
appropriate offsets, since the hstphot pre-processing task coadd is 
unable to apply shifts before combining i mages. Our photo metry 
in this filter is ~0.3 mag brighter than in lli et al.l (2008b) who 
may have combined the exposures without applying offsets, thereby 
clipping the stellar profiles. 

Ground-based adaptive optics observations of SN 2008ax 
where taken on 2008 April 3 using Altair/NIRI on the 8.1-m Gem- 
ini (North) Telescope. These observations were carried out as part 
of programs GN-2008A-Q-28 (PI: Crockett) and GN-2008A-Q-26 
(PI: Gal- Yam). Short exposures (to keep the SN counts linear) to- 
talling 1200 sec on-source were taken in the A^-band using the SN 
as the natural guide star (NGS). The images were reduced, sky- 
subtracted and combined using the NIRI reduction tools within the 
Iraf gemini package. The final reduced image is of very high qual- 
ity showing near-diffraction limited resolution of 0.09" and Strehl 
ratios of up to ~10 percent (Figure 1). The pixel scale of the NIRI 
f/32 camera was 0.022" providing excellent sampling of the PSF. 



3 ANALYSIS AND RESULTS 

A transformation between the pre- and post-explosion coordinate 
frames was derived in order to identify the precise position of 
the SN on the archival HST images. The positions of 35 point 
sources common to both the HST F814W and the Gemini K- 
band frames were input to the Iraf task geomap, which calcu- 
lated a general geometric transformation with a rms error of ±22 
mas. The position of the SN in the Gemini AO image was trans- 
formed to the HST F814W coordinate frame yielding a pixel po- 
sition of [676.08,307.13] on the WF3 chip. Figure 1 shows the 
aligned pre- and post-explosion images, all of which are centred 
on the SN position. A source is visible at the SN site in the 
F450W, F606W and F814W pre-explosion images, although noth- 
ing is seen in the F300W frame (not shown in Figure 1). Its posi- 
tion in the F814W image was measured using 5 different methods, 



DAOPHOT PSF fitting, hstphot and the three centring algorithms 
within DAOPHOT - centroid, gaussian and ofilter. The mean pixel 
position was [676. 1 1,307.17] with an uncertainty of ±12 mas. With 
a measured displacement of 5 mas and a total astrometric uncer- 
tainty of ±25 mas, this source can be considered co incident with SN 
2008ax. This is the same object as suggested by iLi et af (2008a) 
to be the possible progenitor. However, high resolution images are 
necessary to reduce the positional error, and ambiguity, as much as 
possible. 

A slight pixel offset between the F450W and F814W images 
as measured by the Iraf task crosscor is consistent with the offset in 
the position of the proposed progenitor, confirming that this is the 
same source in both filters. The F606W image was taken ~7 years 
earlier with a pointing entirely different from the F450W/F814W 
observations. Positions of 35 point sources common to the F606W 
and F814W images were input to geomap, resulting in a transfor- 
mation with an rms error of ±16 mas. The F814W object position 
was transformed to the F606W coordinate frame and found to be 
coincident within the astrometric uncertainties. 

hstphot flight system magnitudes of the pre-explosion source 
were measured as F450W = 23.66 ± 0.10, F606W = 23.36 ± 0.10, 
F814W = 22.63 ±0.10. From pixel statistics at the SN position we 
estimate a 3<x detection limit of F300W=22.9. 



4 DISTANCE, REDDENING AND METALLICITY 

Estimates of the distance to NGC 4490 and the reddening towards 
the SN progenitor are taken from a compa nion paper on the evolu- 
tion of SN 2008ax |Pastorello et all d2008h . Pastorello et al. derive 
a mean distance modulus of p = 29.92 ± 0.29 and a reddening of 
E(B - V) = 0.3 estimated from the equivalent widt h (EW) of the 
Na I doublet and the relations o f lTurattoet ail [2003 ). We note that 
IChornock etal] d2008l) find a colour-excess of E(B - V) = 0.5 de- 
spite measuring the same EW as Pastorello et al. 

A recent study has estimated the oxygen abundances of 
NGC 4490 from SPSS DR5 spectra of 5 individual Huregions 
( Pilvugin & Thuan 20q3), none of which are close to the SN po- 
sition. Values between 8.3 - 8.5 dex (on the scale of 12 + log O/H) 
are measured around Rq/R2s ~ 0.2, while SN 2008ax is at a galac- 
tocentric distance of 0.3. Pilyugin & Thuan calibrations give abun- 
dances which are slightly, but systematically, lower than previous 
studies. The Hn region abundances measured by these authors at 
the position of SN 2008ax appear similar to those they measure 
in the central regions of M10 1, for which the best recent estimates 
suggest solar-like values (e.g. jBresolin et al.l2007l) . Hence we con- 
clude that the progenitor of SN 2008ax most likely had solar-type 
metallicity and is bracketed from below by LMC-like values (i.e. 
8.3 dex) taking the Pilyugin & Thuan data at face value. 



5 THE PROGENITOR OF SN 2008AX 

These pre-explosion images constitute the most detailed informa- 
tion on a SN progenitor since th e identification of the precursor of 
SN 2003gd dSmartt et aT] 120041) . and the first for a peculiar type 
lib SN since SN 1993 J. In the case of SN 1993 J the progenitor 
was identified as a K-type supergiant with an UV and B-band 
excess interpreted as a nearby hot, mas sive early B-type star 
dAldering et al] 1 19941 : iMaund et al.l |2004|) . Interaction with this 
proposed binary com panion is used to ex plain the stripped nature 
of the SN progenitor. iRvder et alj d2006l) found a point source at 
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Figure 1. Pre- and post-explosion images of the site of SN 2008ax in NGC 4490. Each frame is centred on the SN position and oriented such that North 
is up and East is to the left. The position of SN 2008ax is indicated by the cross hairs, (a) Pre-explosion HST+WFPC2 F450W image, (b) Pre-explosion 
HST+WFPC2 F606W image, (c) Pre-explosion HST+WFPC2 F814W image, (d) Post-explosion Gemini+Altair/NIRI A"-band image in which the SN is 
clearly visible. An object coincident with the SN position is visible in all three pre-explosion frames. The HST+WFPC2 F300W pre-explosion image, in 
which this object is not detected, is not shown here. 



the site of the type lib SN 2001ig ~1000 d post-explosion which 
they identify as the blue supergiant companion to the putative 
progenitor. Given the type lib classification of SN 2008ax one 
might assume a similar binary progenitor system to be appropriate. 
Below we compare the pre-explosion photometry and magnitude 
limits with several progenitor scen arios. Using ext i nction s derived 
with the AJE(B - V) relations of IVan Dvk et alj jl999l) , we find 
absolute photometry for the progenitor of M F 4 S0W = -7.55 ± 0.31, 
M F606W = -7.43 ± 0.31, Mptuw = -7.86 ± 0.31 (Poisson noise = 
±0.10, /j m - = ±0.29) and a 3cr detection limit of Mf3ootv=-8.6. 

Single massive supergiant pr ogenitor. The spectral ene rgy distibu- 
tions of supergiant stars from lDrilling & Lan dolt (2000) were con- 
verted from UBVRI t o HST flight system magnitudes using the 
colour corrections of Maund & Smartt (2005). Comparing these to 
the observerd SED of the pre-explos ion source we f ound no sin- 
gle star solution, in agreement with lLi et al . (2008b). This is not 
surprising given that the source is simultaneously blue in (F450W- 
F606W) and red in (F606W-F814W). Even when applying arbi- 
trary levels of extinction we were unable to fit a single supergiant 
SED to the observed photometry. 

Superposition of two supergiant s tars. That SN 2008ax transi- 
tioned from a type II to a type lb SN jPastorello et al.l2008h implies 
that the progenitor had lost all but a small fraction (< few 10~'M G ) 
of its H-rich envelope. Previous studies suggest that the type lib 
SNe 1 993J and 2001ig arose from massive binary progenitor sys - 
tems jAldering et ai]|l994l ; lMaund et alj[2004l : iRvder et all 2006), 
the binary orbit in each case being close enough to allow interac- 
tion. For SN 2008ax we do not restrict ourselves to an interacting 
binary scenario. The PSF of the proposed progenitor is consistent 
with that of a single star; however, its FWHM is some 0. 15", which 
translates to ~ 6pc at the distance of NGC 4490. The source might 
actually be two stars separated by as much as 3 pc and such widely 
separated objects would not interact, evolving instead as two sin- 
gle stars. The progenitor in this case would have to lose almost its 
entire hydrogen envelope through strong stellar winds prior to ex- 
ploding, a feat which theo ry suggests is only poss ible for stars more 
massive than -25-30 M Q feldridge & Toutll2004h . 

Pairs of supergiant SEDs from Drilling & Landolt (2000), 
converted to the HST VEGAMAG system (see above), were fit- 
ted to the source photometry with the same value of extinction 
being assumed for both components. The best fit was a B8 su- 



pergiant of log(LIL a ) »5.15 combined with a M4 supergiant of 
log(L/L Q ) »4.20. This implies the initial mass of the more evolved 
M-type supergiant (~ 8 - 10M o ) is much less than that of its B-type 
companion (~ 25M G ). Assuming the stars are coeval one would ex- 
pect the more massive object to be the most evolved. Even if one 
argues that the more massive object is evolved, the lower mass star 
should not be observed as a red supergiant as its main sequence 
lifetime is ~4-5 times longer. We also attem pted to fit a Wolf -Rayet 
(WR) + supergiant interacting binary (e.g. jRvder et al. 2004) using 
model WR colours (see discussion of single Wolf-Rayet progenitor 
for details) and found a similar result; the WR star was around 50 
times more luminous than its evolved red/yellow supergiant com- 
panion. Neither of these two-star models is self-consistent. 

The extinction towards an assumed companion star was var- 
ied, while retaining E(B - V) = 0.3 for the progenitor. Since this 
implies that the stars are separated by some considerable distance, 
we cannot invoke binary interaction. Arbitrarily reducing the com- 
panion extinction to zero we found a range of fits where the com- 
ponents had quite similar luminosities/masses, the progenitor being 
of later spectral type. The relative evolutionary timescales are not 
inconsistent in this case; however the implied main sequence mass 
of 10-14 M Q is much too low for the progenitor to have lost most 
of its hydrogen envelope through wind-driven mass-loss. 

It is possible that several other stars migh t lie within the PSF as 
was the case for the progenitor of SN 1993 J jAldering et al]| 19941 ; 
Maund ct "all l2004l) . Removing the flux contributed by these stars 
might allow us to fit a consistent interacting binary model. 

Young stellar cluster. Alternatively the source might be an un- 
resolved, young stellar cluster in which the SN progenitor was 
embedded, provided the cluster diameter is less than ~6 pc (PSF 
FWHM). This is not unreasonable as compact stellar cl usters ob- 
served in M51 have median e ffective radii of 2-4 pc jLee et all 
|2005| ; IScheepmaker et al.ll2007D . Previous studies have attempted 
to place constraints on SN progenitors by estimating the ages and 
main se quence turn-off masses of t heir host clus ters (e.g., SN 
2004d j - iMafz-Apellaniz et alj[2004al; SN 2007gr -ICrockett et al.1 
|2008| ; GRB 030329/SN 2003dh - lOstlin et al.|[2008l) . Here we at- 
tempt a similar analysis ass uming that the pre- explosion source 
is a compact, coeval cluster. iBastian et al. I d2005h suggest sources 
brighter than M v =-&.6 are most likely star clusters. Our source 
is somewhat less luminous (M F606W = -7.43) so cannot imme- 
diately be characterised as such. We used CHORIZOS (version 
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2.1.4) <Mafz-Apellanij 2004b). 

3.x 2 minimisation code, to fit Star- 



burst ^ (Leither er et al.M 1999T) stellar populations to our observed 
photometry. The model SEDs were of solar metallicity and as- 
sumed a Salpeter initial mass function with an upper mass cut- 
off of 1OO M . CHORIZOS solved for clus ter age and extinc- 
tion using a lCardelli. Clayton & Mathis] dl989h extinction law with 
R v =3.l. The best fit was a cluster with an age of ~20 Myr (turn- 
off mass ~ 12M G ), a total mass of ~ 1800M Q and extinction of 
E(B - V) = 0. 15, although a range of almost equally significant so- 
lutions was found with ages of 8-28 Myr and E(B - V) = 0.0 - 0.2. 
This degeneracy was due to the lack of information at UV and 
NIR wavelengths. Since none of these 'best fits' had an extinc- 
tion consistent with our measured value, we forced CHORIZOS 
to adopt E(B - V) = 0.3. The solution in this case was a 42 Myr 
cluster (turn-off mass ~ 8M Q ); however, the model SED was under- 
luminous in F450W and over-luminous in F606W compared with 
the observed photometry. 

Given the relatively low mass estimates for these clusters, it is 
perhaps not surprising that we do not find a satisfactory match be- 
tween the models and our observed SED. The model SEDs assume 
an infinite number of stars in the cluster and therefore a fully pop- 
ulated IMF; a reasonable approximation for more massive clusters. 
For low mass clusters, where the IMF is more sparsely populated 
and the total luminosity is much lower, even a single massive star 
can significantly affect the observed SED; and the stochastic nature 
of the masses of the stars, especially the most massive stars, can 
result in a large dispersio n in the integrated photometry of clus- 
ters o f similar masses (e.g.. lCervino & Luridianal2004l ; |jamet et al.l 
l2004h . In an attempt to mimic these effects we created a grid of 
'cluster+massive star' models by combining the flux from a sin- 
gle, evolved star (of mass 10-30 M Q and spectral type WR/09-M5) 
with coeval (~6 to 27 Myr) Starburst99 stellar populations. The 
single massive star was assumed to be the SN progenitor. These 
hybrid models were fitted to the observed photometry adopting a 
value of E(B - V) = 0.3. The resulting fits yielded cluster masses 
of ~ 1500 - 3000M o , main-sequence masses for the progenitor 
of 10-14 M Q , and colours consistent with OB-type stars and WR 
objects that have retained a fraction of their H envelope. Stars of 
such low initial mass could not have produced a H-poor progeni- 
tor through wind-driven mass-loss. Hence we conclude that, if em- 
bedded in a cluster, the progenitor star must have been stripped in 
an interacting binary system. Interestingly these initial masses are 
similar to that estimated for the progenitor of SN 1993J (~15 M Q ), 
although it exploded as a cooler, K-type star. 
Single Wolf-Rayet star. A progenitor star of sufficient initial mass 
could have lost most of its H-rich envelope through radiatively 
driven mass loss. WR stars, exposed He cores of massive stars, 
can form in this way. Helium-rich WN stars have been proposed 
as possible progenitors of type lb SNe and H-rich supe rgiants have 
been shown to be the progenitors of type II SNe (e.g.. lSmartt et al.l 
l2004l;lLietal]|2006h . A progenitor intermediate of these types of 
star would be consistent with the strange evolution of a type lib 
SN. 

We compared our pre-explosion photometry with solar metal- 
licity models of single stars of initial masses between 20 and 32 
M Q . Our model c olours and absolute magnitu des were obtained us- 
ing the method of lLeieune & Sch aerei (2001) to perform synthetic 
photometry on the Ca mbridge STARS stellar evolution models 
dEldridge & Vink!l2006l). This process is described in more detail 
by Eldridge, Mattila & Smartt (2007). although here wc have also 
calculated s ynthetic photometry for WR stars using the model WR 
spectra of Grafener et al. (2002) and Hamann & Grafener (2003). 



Photometry was calculated using both the Johnson-Cousins and 
HST WFPC2 photometric systems, so that it was possible to com- 
pare the models directly with BVI photometry of observed WR 
stars and WFPC2 photometry of the proposed progenitor respec- 
tively. Figure 2 shows BVI colour-magnitude and colour-colour 
plots of these models along with WN and WNL stars in M3 1 
and the LMC dMassev et al] 120051 ; iMassevI l2002h . The photom- 
etry of the proposed progenitor, transf ormed to BVI magnitudes 
(Holtz manetal]|l995l ; Dolphin 2000b), is plotted with a solid 
square and its colours fit well with those of the observed WR stars. 
We note that the observed WR stars show a considerable dispersion 
in colour, particularly in V-I, when compared to the model tracks. 
The strong mass-loss and emission-line dominated spectra of these 
objects lead to an in herent spread in their broad-band colours (e.g., 
Ham ann etal]|l995h . Given the ambiguous nature of the object's 
SED our transformation WFPC2 to Johnson-Cousins magnitudes 
maybe somewhat uncertain, but comparison of the flight-system 
magnitudes with the WFPC2 model photometry shows a very simi- 
lar relationship between the pre-explosion source and the models to 
that shown in Figure 2. The progenitor lies beween the endpoints of 
the 27 and 28 M Q tracks. Remarkably these model stars end their 
lives while transitioning between the supergiant and WR phases, 
precisely the type of progenitor one would expect to produce a type 
lib SN. Final hydrogen masses are between ~ IP - 1 and 10~ 3 M H . 
These would be classed as WNL or WNH stars dSmith & Contil 
120081) . 

A potential pitfall for our single WNL/WNH progenitor is 
the mass of its carbon/oxygen (C/O) core. The final stellar mass 
is 11-12 M , with a C/O core mass of ~9 Mb; a factor of tw o 
higher than for the pr ogenitor of SN 1993J dMaund et al]|2004l) . 
IPastorello et alj (2008) show that the bolometric light curves and 
ejecta velocities of SNe 2008ax and 1993J are very similar, in- 
ferring that ejecta masses, and hence progenitor masses were also 
similar. They suggest the more massive C/O core arising from our 
single star model would have produced a broader light curve. If 
these two progenitors actually had core masses and ejecta masses 
different by a factor of two, then the fact that they show similar 
evolution would have important implications for modelling of SN 
lightcurves. Alternatively, the ejecta masses could be similar if we 
assume our WNL/WNH progenitor produced a more massive com- 
pact remnant (~5 M Q ) than was the case for SN 1993J. The rem- 
nant formed would be a black hole, in agreeme nt with the predic- 
tions of stellar models in this initial mass range dHeger et al.ll2003l; 
lEldridge & Toutl2004h . 



6 CONCLUSIONS 

We have determined the position of SN 2008ax to within +22 mas 
in pre-explosion HST observations of NGC 4490 and identified a 
coincident source in three filters. It is unclear whether this object is 
the SN progenitor alone, or a blend of it and several other stars. We 
conclude that it is not a single supergiant star, but that an OB/WR 
star progenitor (of main-sequence mass 10-14 M Q ) in an interact- 
ing binary is possible so long as we include a significant flux con- 
tribution from an unresolved stellar cluster. We also find a single 
star progenitor model consistent with the pre-explosion photome- 
try, where an initially very massive star (~28 M Q ) loses most of 
its H-rich envelope prior to exploding as an 11-12 M G WNH star. 
However, the light curve of SN 2008ax would suggest a lower pro- 
genitor mass ( Pasto relloetai] |2008). In several years time, when 
the SN has faded, further high resolution imaging may help us bet- 
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Figure 2. Comparison of the pre-explosion source photometry with STARS stellar evolution models I Eldri dge & Vinki 2006). The blue lines are non-rotating 
solar metallicity models between 20-32M Q . The thick black line is the 27M Q model; the dashed line the 28M Q model. The red, green and bl ue squares indicate 
the endpoints of the 27, 28 and 29M Q models respectively. The black asterisks are observed WN and WNL stars in M31 and the LMC iMassev et al. 2005; 
Massey 2002), which have been corrected for Galactic extinction. The possible SN progenitor is marked by the black square, with the hollow square its position 
before correcting for extinction. 



ter determine the nature of the progenitor star. If a source is still 
visible, subtraction from the pre-explosion data could help to iso- 
late the progenitor SED. If the source has disappeared, it could have 
serious implications for the interpretation of ejecta masses from SN 
lightcurve and spectral modelling. 
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